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We describe a novel scheme for low-noise phase-insensitive linear amplification at microwave fre- 
quencies based on the Superconducting Low-inductance Undulatory Galvanometer (SLUG). Direct 
integration of the junction equations of motion provides access to the full scattering matrix of the 
SLUG. We discuss the optimization of SLUG amplifiers and calculate amplifier gain and noise tem- 
perature in both the thermal and quantum regimes. Loading of the SLUG element by the finite 
input admittance is taken into account, and strategies for decoupling the SLUG from the higher- 
order modes of the input circuit are discussed. The microwave SLUG amplifier is expected to achieve 
noise performance approaching the standard quantum limit in the frequency range from 5-10 GHz, 
with gain around 15 dB for a single-stage device and instantaneous bandwidths of order 1 GHz. 



I. INTRODUCTION 

The rapid development of superconducting quantum 
electronics has motivated a search for near quantum- 
hmited microwave amplifiers for the low-noise readout of 
qubits and linear cavity resonators. It was long ago rec- 
ognized that the dc Superconducting QUantum Interfer- 
ence Device (dc SQUID) can achieve noise performance 
approaching the fundamental quantum hmit imposed on 
phase-insensitive linear amplifiers: namely, the amplifier 
must add at least half a quantum of noise to the sig- 
nal it amplifies [1]. Yet while the SQUID is in principle 
capable of amplifying signals at frequencies approaching 
the Josephson frequency (typically of order tens of GHz), 
it remains challenging to embed the SQUID in a 50 fi 
environment and to provide for efficient coupling of a 
microwave signal to the device. Recently, it was shown 
that near quantum-limited performance can be achieved 
with a microstrip SQUID amplifier, where the input coil 
is configured as a microstrip resonator with the SQUID 
washer acting as a groundplane [2]. The noise tempera- 
ture of a microstrip SQUID amplifier cooled to millikelvin 
temperatures has been measured to be 47 ± 10 mK and 
48 ± 5 mK at frequencies of 519 MHz and 612 MHz, re- 
spectively, more than an order of magnitude lower than 
the best semiconductor amplifiers available and within a 
factor of 2 of the quantum limit [3l |4] . However, efforts 
to extend the operating frequencies of these amplifiers 
into the gigahertz range are hampered by the fact that 
reduction of the length of the input resonator is coupled 
to reduction of the mutual inductance between the input 
coil and the SQUID [5|. Alternative approaches have in- 
cluded the integration of a high-gain SQUID gradiometer 
into a coplanar waveguide resonator at a current antinode 
017]. 

The current study was motivated by the development 
of a new device configuration that enables the efficient 
coupling of a GHz-frequency signal to a low-inductance, 
high gain SQUID that should achieve noise performance 
approaching the standard quantum limit. The gain ele- 
ment is more properly termed a SLUG (Superconducting 



Low-inductance Undulatory Galvanometer), as the sig- 
nal is not coupled to the device inductively, but rather 
injected directly into the device loop as a current [8]. 
The low-inductance design is straightforward to model 
at microwave frequencies, and the SLUG is readily in- 
corporated into a microstrip line in such a way that the 
modes of the SLUG element and the input resonator re- 
main cleanly resolved, greatly simplifying analysis of the 
circuit. In what follows we present a comprehensive the- 
oretical study of the gain and noise performance of the 
SLUG microwave amplifier. Our goals are to clearly spell 
out the design tradeoffs, to outline a clear path to device 
optimization, and to identify the fundamental limits to 
performance. 

As we shall see, the scattering parameters of the SLUG 
are very similar to those of the more familiar symmetric 
dc SQUID, apart from a trivial shift in flux bias that 
arises from the asymmetric division of bias current be- 
tween the two arms of the SLUG. However, while it is 
straightforward to fabricate a low-inductance (~ 10 pH) 
SLUG and to embed the device in a 50 fi environ- 
ment, it is challenging to engineer a clean, purely in- 
ductive coupling to a conventional dc SQUID at mi- 
crowave frequencies. For this reason we have chosen 
to focus our discussion of microwave amplifiers on the 
SLUG geometry. In this manuscript we will not con- 
sider phase-sensitive amplifiers based on parametrically 
modulated Josephson junctions operated in the super- 
current state [9j Ho] . There has been significant recent 
development of low-noise Josephson parametric ampli- 
fiers [TTHI3]j including such milestones as squeezing of 
vacuum noise '14] and observation of quantum jumps 
in a superconducting qubit J15j . Because these ampli- 
fiers squeeze the input state, they can achieve added 
noise numbers for one field quadrature that are below 
the standard quantum limit JTHIIIT]. Moreover, these de- 
vices operate with negligible dissipation, circumventing 
practical problems associated with hot-electron effects 
that are intrinsic to devices that operate in the finite- 
voltage regime. In related work, there have been efforts 
to develop low-noise phase-insensitive amplifiers based 



on parametrically modulated junctions in a ring modu- 
lator configuration [18]. Broadly speaking, advantages 
of the Josephson parametric amplifiers include unsur- 
passed noise performance and ease of fabrication, while 
potential disadvantages relative to SQUID-based dissipa- 
tive amplifiers include modest gain-bandwidth product, 
limited dynamic range, and increased complexity of op- 
eration. Ultimately we suspect that there is a place in 
the superconducting quantum optician's toolbox for both 
ultralow noise phase-sensitive parametric amplifiers and 
robust, broadband phase-insensitive amplifiers operating 
near the standard quantum limit. 

This paper is organized as follows. In Section II we 
introduce the circuit models of the symmetric dc SQUID 
and the SLUG. In Section III we calculate the dc char- 
acteristics of the devices. In Section IV we evaluate 
SLUG scattering parameters, and examine the maximum 
achievable gain over the range of device parameters. Sec- 
tions V and VI present an analysis of noise properties in 
the thermal and quantum regimes, respectively. In Sec- 
tion VII we describe the design and performance of prac- 
tical SLUG amplifiers for GHz frequency operation, and 
in Section VIII we discuss amplifier dynamic range. In 
Section IX we describe the effect of the finite admittance 
of the input circuit on device characteristics, gain, and 
noise, and in Section X we discuss hot-electron effects. 
In Section XI we present our concluding remarks. 



II. DEVICE MODEL 



To make contact with the earlier numerical studies of 
Tesche and Clarke ^9,, we begin by considering the fa- 
miliar symmetric dc SQUID, shown in Fig. [Tk. The gain 
element consists of two overdamped Josephson junctions 
embedded in a superconducting loop with inductance L. 
The junctions (with gauge invariant phases <5i^2) have 
equal critical currents /q, self-capacitances C, and shunt 
resistances R. The superconducting loop is formed from 
two equal branches with inductance i/2; we neglect the 
mutual inductance between the branches. A dc bias cur- 
rent /{, and bias flux $;, establish a quasistatic operating 
point, and the signal is injected into an input coil that is 
coupled to the SQUID loop with mutual inductance M . 
The currents through the junctions are given by 



h = Iq sm (5i H — ^ ^d -JT 

R at 

T r ■ X (^2 - K.2) , ^ dV2 ,^, 



where Vn.i and Vn,2 are noise voltages associated with the 
resistive shunts, and where the voltages Vi_2 are related 
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FIG. 1. Device geometries, (a) Symmetric dc SQUID, (b) 
Symmetric SLUG. 



to the junction phases by the ac Josephson relation: 

<i>o dSi 



U2- 



27r dt 

<^Qd62 



(2) 



Here, $o = h/2e is the magnetic flux quantmn. The 
SQUID loop supports a circulating current J given by 



J = 



h-h 



(3) 



The voltage across the device is given by 

V 



^ 2 dt 



V2 + 1, 



L dh 



2 dt 



(4) 



The circulating current and the junction phases are re- 
lated to the total flux in the loop ^t as follows: 



~ 27r 



+ LJ 
(62-61). 



(5) 



We introduce dimensionless variables i, v, (/>, and 9, 
defined as follows: i = I/Iq, v = V/IqR, <j> = ^/^o, 
and 6 = t / [<I>o/(27r/oi?)]. In addition, we introduce the 
dimensionless reduced inductance /3l = 2/oL/$o and the 
damping parameter f3c = {2ti/^q)IqR?C. The equations 
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FIG. 2. I-V characteristics of (a) symmetric dc SQUID and 
(b) symmetric SLUG for $6 = 0<I>o (black), 0.25 $0 (blue), 
and 0.5 <l?o (red). The device parameters are (3l ~ 1 and 
Pc = 0.8. 



FIG. 3. V-$ characteristics of (a) symmetric dc SQUID and 
(b) symmetric SLUG for It = 1.8 Jo (black), 1.9 7o (blue), and 
2.0/0 (red). The device parameters are /3l ~ 1 and Pc ~ 0.8. 



of motion for the junction phases are written as 
ib , 62- Si ~ 2TT(f)b 
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sin(52 - (^2 +'y„,2- (6) 



The quasistatic output voltage and circulating current 
are given by 
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Vout 
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In the SLUG geometry of Fig. [T]d, the device loop 
is formed from two superconducting traces separated by 
a thin dielectric layer, and the input signal is injected 
directly into one of the traces. In the case where the 
SLUG is integrated into a microstrip transmission line, 
the device is realized in three metallization steps (corre- 
sponding to the circuit groundplane and the two arms of 
the SLUG), with two dielectric thin films separating the 
metal layers. The mutual inductance between the arms 
of the device is of order the self-inductance of the arms, 
and must be taken into account. Below for the sake of 
simplicity we consider the case where the two dielectric 
films separating the superconducting layers are of equal 



thickness, resulting in branch inductances 2L and L with 
mutual inductance L. The total inductance of the device 
is then L, and the mutual coupling M of the input current 
/$ to the device loop is also L. We refer to this config- 
uration as the symmetric SLUG. The total flux through 
the device becomes 






(9) 



We write the dimensionless equations of motion for 61^2 
as follows: 

Pcdi = ^^ «0 - smdi - di + vn,i 



Pc52 



(52 — <5i — 27r0(, 



+ jfc + «0 - sin (52 - §2 + vn^2- 

(10) 



The output voltage and circulating current are given by 



Vout 



(11) 



J = -^ (<52 - <5i - 27r<^b) - z^/2. (12) 

TrpL 

To operate the SQUID or the SLUG as an amplifier, one 
chooses lb and <&;, to establish a quasistatic operating 
point where the transfer function V$ = dV/d^ is large. 
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the resulting device is immune from stray reactances and 
straightforward to model at microwave frequencies. The 
required critical current density is 10 kA/cm^, within the 
reach of standard Nb-AlO^j-Nb technology. While Joule 
heating in the shunt resistors is significant, the addition 
of large-volume normal-metal cooling fins should allow 
equilibration of the shunt resistors at temperatures below 
100 mK (see Section X). Lower device inductances would 
require uncomfortably high junction critical currents to 
achieve comparable device performance, and fabrication 
yield and Joule heating of the shunts would become prob- 
lematic. On the other hand, a significantly larger SLUG 
inductance would provide less gain and complicate the 
microwave engineering, owing to the larger device dimen- 
sions. 



IV. SCATTERING PARAMETERS 



FIG. 4. Forward transfer function V<s> of SLUG circuit versus 
quasistatic bias flux for bias currents lb = 1.8 7o (black), 1.9 Jo 
(blue), and 2.0 Jq (red). The device parameters are /3l = 1 
and /3c = 0.8. 



In both cases, the device acts as a transimpedance el- 
ement: the input signal is coupled to the device as a 
current, and the output signal is coupled from the device 
as a voltage. 



In order to optimize SLUG amplifier design, it is nec- 
essary to understand the forward transfer function and 
the complex input and output impedances of the de- 
vice. To extract these from our model, we apply ap- 
propriate test signals and probe the complex response 
at the excitation frequency, chosen to be a small frac- 
tion of the Josephson frequency loj /2tt. The forward 
transimpedance Vi = dV/dl is readily derived from the 
SLUG flux-to- voltage transfer function V$: 



Vi =MU$, 



(13) 



III. DC CHARACTERISTICS 



Equations ([6]) and ( 10 ) were numerically integrated us- 
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time 



ing a 4th order Runge-Kutta solver for N 
steps A9 over a range of bias points. In Figs. [2^-b we 
show the I-V characteristics of the symmetric dc SQUID 
and the symmetric SLUG with /3i = 1 and (3c — 0.8; in 
Figs. [3^-b we show the V-^ characteristics of the same 
devices. For bias near 1.9 /q, the peak-to-peak voltage 
modulation is around 0.5 IqR. We observe that the dc 
characteristics of the SLUG closely match those of the 
symmetric dc SQUID, apart from a shift in flux bias point 
that arises from the asymmetric division of the SLUG 
bias current. Similarly, we have found that the scatter- 
ing parameters and noise properties of the SLUG and the 
SQUID are closely matched, apart from this bias shift. 
Therefore for the sake of simplicity we choose to focus in 
the following on the device characteristics of the SLUG 
alone. We will consider the following set of SLUG param- 
eters: ^i = 1, /3c = 0.8, L = 10 pH, and C = 50 fF, cor- 
responding to a junction with critical current 100 /iA and 
area around 1 /zm^. Several considerations lead us to this 
choice. First, inductances of order 10 pH are realized in 
a reliable, controlled way using the SLUG geometry, and 



where again we have M = L for the case of the symmetric 
SLUG. In Fig. |4] we plot SLUG V$ versus flux over a 
range of current bias points for /3i^ ~ I and f3c = 0.8. 

Next we consider the input return loss. The SLUG 
input is an inductive short to ground at low frequencies, 
and the complex input impedance Zi is frequency depen- 
dent. The input impedance is readily derived from the 
dynamic impedance Z, defined in terms of the flux-to- 
current transfer function J$ = dJ/d^ as follows: 



— Jih ^ -— — 



1 
C 






(14) 



where following Hilbert et al. we have introduced the 
frequency-independent dynamic resistance TZ and dy- 
namic inductance £ [20 . In Figs. [5] and [6] we plot R/TZ 
and L/ C, respectively, for a SLUG with j3i^ ~ 1 and 
Pc — 0.8 over a range of bias points. 

Finally, in Fig. [7]we show the device output impedance 
Ro over a range of bias points. The output impedance 
is real and frequency independent, and the magnitude of 
Ro is of order the junction shunt resistance R. 

For the following discussion, it is convenient to work 
in terms of the bias-dependent dimensionless impedance 
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FIG. 5. R/n versus flux for a SLUG with /3l = 1 and /3c = FIG. 7. SLUG output resistance Ro versus flux for bias cur- 
0.8, for bias currents 7^ = 1.8 Jo (black), 1.9 Jo (blue), and rents h = 1.8 /o (black), 1.9 7o (blue), and 2.0 Jo (red). The 
2.0 7o (red). device parameters are Pl ~ 1 and /3c ~ 0.8. 
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FIG. 6. L/C versus flux for a SLUG with Pl = 1 and /3c = 
0.8, for bias currents h = 1.8 Jo (black), 1.9 Jo (blue), and 
2.0 7o (red). 
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FIG. 8. Maximum achievable power gain Gm for a SLUG 
amphfier versus flux for bias currents It = 1.8 7o (black), 
1.9 7o (blue), and 2.0 7o (red). The device parameters are 
/3l = 1, /3c = 0.8, 7/ = 10 pH, and C = 50 fF; the operating 
frequency is 5 GHz. 



parameters pi^o, defined as follows: 

(wM)2 



Ri — Pi 



R 



Ro = PoR (15) 

From the definition of TZ it follows that pi — R/TZ. As we 



will see, amplifier gain, bandwidth, and noise properties 
depend sensitively on pi and po- 

Power gain of the device is maximized when appropri- 
ate conjugate matching networks are employed to couple 
the signal to and from the device. The maximum avail- 



able power gain Gm is given as follows: 



G,n — 



I?R, 



(16) 



where li is the input current and Vo is the output voltage. 
Using Eq. (15 1, we find 



^rn 



'ipiPo 



w J 



(17) 



In Fig. [8] we plot Gm for the symmetric SLUG with 
/3i = 1, /3c = 0.8, L = 10 pH, and C = 50 IF for an op- 
erating frequency of 5 GHz. Over a broad range of bias 
parameters gain in excess of 20 dB is readily achievable. 
It is important to note, however, that a conjugate match 
to a 50 $7 source does not yield best amplifier noise per- 
formance, due to the mismatch between the real part of 
the SLUG input impedance Ri and the optimal noise- 
matched source impedance, which can be significantly 
larger than Ri. Amplifier optimization therefore involves 
a tradeoff between gain and noise performance, as dis- 
cussed in detail below. 

The bandwidth of the SLUG amplifier will be deter- 
mined by the coupling to the low- impedance input port, 
as the device output is reasonably well-matched to typi- 
cal transmission line impedances. To get a rough idea of 
amplifier bandwidth we consider a 50 SI source impedance 
and assume that conjugate matching at the device input 
is accomplished via a simple quarter-wave transmission 
line section; for simplicity we neglect the imaginary part 
of the SLUG input impedance. The amplifier quality 
factor Q is given by 



(18) 



The bandwidth of an amplifier designed at an operat- 
ing frequency uj/2t: is then uj/2ttQ. For an operating 
frequency around 5 GHz, we find that R is of order 
0.1 f2. Therefore we expect Q of order 10, and ampli- 
fier bandwidths of order hundreds of MHz. For current 
bias lb < 2/o and for a narrow range of fluxes correspond- 
ing to bias points near the supercurrent branch, we find 
that it is possible to achieve extremely high power gain 
(see Fig. Isl). However, the high gains achieved at these 
bias points are due largely to vanishing R; an amplifier 
designed to operate in this regime would have a rather 
small bandwidth. It is important to note that Eq. ( 18 ) 




presents only a rough guideline for the bandwidth rather 
than a fundamental limit. In particular, it is possible to 
obtain a larger bandwidth with no degradation in gain 
by employing either a tapered transmission line matching 
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FIG. 9. Dimensionless SLUG noises (a) 7v, (b) 7.7, and (c) 
yvj versus flux for bias currents h — 1.8 Jo (black), 1.9 /o 
(blue), and 2.0 /o (red). The SLUG parameters are /3l — 1, 
/3c = 0.8, and r = 4 X 10"^ 



section or a multisection input transformer with stepped 
transmission line impedances. We postpone a more de- 
tailed discussion of amplifier bandwidth to Section VII. 




FIG. 10. Circuit for noise analysis. 
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FIG. 11. Real part of the optimal source impedance Rs.opt 
versus flux for bias currents /;, = 1.8 7o (black), 1.9 Jo (blue), 
and 2.0/0 (red). The SLUG parameters are /3l = 1, /3c = 0.8, 
and r = 4 X 10^"". The operating frequency is 5 GHz. 



Following Clarke et al., we introduce the dimensionless 
noise parameters 7y, 7j, and 7yj, such that the voltage 
noise spectral density at the device output is given by 
Sv — 2jykBTR, the circulating current noise spectral 
density is Sj = 2jjkBT/R, and the cross noise spec- 
tral density is Syj — 2jvjkBT; here T is the electron 
temperature of the shunt resistors [5T1 1^ . These noises 



are calculated by solving the Langevin equations (10 1 



The noise spectrum consists of a series of peaks at the 
Josephson frequency and its harmonics; the dimension- 
less noises 7 are evaluated at low frequency / <^ ujj/2it 
where the spectrum is white. The noises 7 do depend on 
the noise parameter T, due to the possibility of satura- 
tion and smearing of the device characteristics at elevated 
temperature. In Fig. [9] we plot the dimensionless noises 
over a range of bias parameters of the symmetric SLUG 
for Pl ^ 1, Pc = 0.8, and r = 4 X 10"^ this choice 
corresponds to a temperature of 100 niK for a junction 
critical current of 100 /lA. We note that at high bias cur- 
rent /(, ^ -^0: "fv,J approach the expected Johnson noise 
limit of 1 for the two shunt resistors in parallel. 

The device noise temperature T„ can be evaluated from 
the circuit shown in Fig. [TO] We assume a noiseless 
source impedance Zg ~ Rg + jX^ and equate the total 
noise of the amplifier to the noise contribution from a 
source resistance Rs at an effective temperature T„. We 
refer all noises to the device output. We find 



AksTnRs „2 

Rf 



2'yvkBTR 



V|M2 



X^ 



2-fjkBT uj^ViM^ 



R 



R't 



X? 



ijvjkBT 



jV^JVPXt 

W 
(19) 



R^ 



Here Rt — Rs + Ri {Xt — Xg + Xi) is the sum of the 
real (imaginary) parts of the source impedance and the 
device input impedance. The noise temperature is thus 
given by 



V. NOISE PROPERTIES IN THE THERMAL 

REGIME T„ 

The Johnson noise of the shunt resistors gives rise to 
a voltage noise at the device output and to a circulat- 
ing current noise in the device loop; moreover, these 
noises are partially correlated, since the circulating cur- 
rent noise couples a flux noise to the loop, which in turn 
yields a voltage noise across the device. To incorporate 
noise in our model, we used a pseudorandom number 
generator to create a gaussian-distributed set of voltages 
vn^i and V]\[^2 with zero mean and variance 2T/A0, where 
we have introduced the dimensionless noise parameter 
r — 27rfcBT//o$o; this choice corresponds to the usual 
white power spectral density S^ = AT for Johnson noise 
in the thermal limit. The simulations were averaged over Rs,opt 

many (^ 100) realizations of the random noise voltages. 



jv{Rl + X?)R 
2 ViAPRs 



2 RRs 



LoXt 



IVJ 



V^Rs 



T. 

(20) 



We use the condition dTn/dXf = to solve for the 
imaginary part of the optimal source impedance. We 
find 
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iv 
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(21) 



Similarly, the condition dTn/dRg = yields the real part 
of the optimal source impedance. We have 



1 fy'S'V f 2 N 
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FIG. 12. Optimal SLUG noise temperature versus flux for 
bias currents Ji = 1.8 Jo (black), 1.9 /o (blue), and 2.0 Jo (red). 
The SLUG parameters are /3i = 1, /3c = 0.8, and L = 4 x 



10" 



The operating frequency is 5 GHz. 



VI. NOISE PROPERTIES IN THE QUANTUM 
REGIME 

At sufficiently low temperature, the zero-point fluctu- 
ations of tlie resistive shunts are expected to make the 
dominant noise contribution. The full expression for the 
spectral density of voltage noise produced by the resis- 
tors is written as 2hfRcoth(hf/2kBT). We have cal- 
culated the added noise of the symmetric SLUG in the 
zero-temperature limit, where the voltage spectral den- 
sity of the shunt resistors becomes 2hfR. We generate 
a single-sided quantum spectral density by digitally fil- 
tering gaussian white noise. Using the quantum noise 



as a driving term in the Langevin equations (10), we 



evaluate the voltage power spectral density Sy (/) at the 
device output, the circulating current spectral density 
S,j{f), and the cross spectral density Svj{f); in Fig. 
[13] we plot these noises versus flux for various bias cur- 
rents. Once again, the device noise temperature T„ can 
be evaluated from the circuit of Fig. [TOJ We assume 
a zero-temperature source impedance Z^ = Rs + jX^ , 
and equate the total noise of the amplifier to the noise 
contribution from a source resistance it!., at a finite effec- 
tive temperature T„. The amplifier noise temperature is 
obtained from the relation 



For bias points where V$ is highest, we have the following 
approximate expression for Rs,opt'- 



R 



s,opt 



1 V^ 

IvPi w 

IvR 



[IvlJ 



I^jY^R. 



{ivlj - Ivj) 



1/2 



(23) 
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(24) 



In Fig. 12 we show the optimal noise temperature Tn,opt 
for a SLUG amplifier over a range of bias points at an op- 
erating frequency uj/2t: = 5 GHz. Note that every point 
in these plots corresponds to a different realization of the 
input matching network; in Section VII we will examine 
the bias- and frequency-dependent noise temperature of 
SLUG amplifiers operated with a fixed input network. 
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(25) 



In Fig. 11 we plot Rs,opt/Ri versus flux for various 

bias currents. For typical device parameters, we have 

Rs.opt S> Ri- For this reason, it is not possible to achieve 

a simultaneous power match and noise match. It is go that 

worthwhile to note, however, that the ratio Rs.opt/Ri 

scales with frequency as lo~^, facilitating simultaneous 

attainment of high gain and good noise performance at 

higher operating frequencies. 

When the signal is coupled to the device via a source 
with optimal impedance Rs,opt -^ jXg^opt, the amplifier 
noise temperature becomes 



Alternatively, we can express the noise contribution of 
the device in terms of an added number of noise photons 
n, where n and T„ are related as follows: 



coth {hf/2kBTn) = 2n + l, 



(26) 



2hfRs 



Sv Rt + X^ 
^ U|Af2 
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(27) 



The optimal source impedance Zg^opt = Rs,opt + jXs^opt 
is obtained from the relations dn/dXt = and 
dn/dRs = 0. The imaginary part of the optimal source 
impedance is given as follows: 
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ujV^M^ - X, 



(28) 



Similarly, the real part of the optimal source impedance 
is written 
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FIG. 14. Real part of the optimal source impedance Rs,opt in 
the quantum regime versus flux for bias currents h — 1.8 Jo 
(black), 1.9 Jo (blue), and 2.0 Jo (red). The operating fre- 
quency is 5 GHz and the SLUG parameters are Pl = 1, 
Pc = 0.8, L = 10 pH, and C = 50 fF. 



In the limit V$ 3> w, we find 
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(30) 



In Fig. 14 we plot Rs,opt/Ri in the quantum regime ver- 



sus flux for a range of bias currents. 

For the optimally matched source, the added number 
of noise photons is given by 
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(31) 



In Fig. 15 we plot Uopt versus flux, for various current 
biases. We see that for an appropriately noise-matched 
source the SLUG approaches a noise level that is close to 
the standard quantum limit usql = 1/2, the minimum 
achievable added noise for a phase-insensitive linear am- 
plifier [T6] . 



FIG. 13. Quantum noises (a) Sv, (b) Sj, and (c) Svj versus 
flux for bias currents /j, — 1.8 Iq (black), 1.9 7o (blue), and 
2.0 7o (red). The SLUG parameters are Pl = I, Pc = 0.8, 
i = 10 pH, and C = 50 fF. 



VII. AMPLIFIER DESIGN 

The above analysis demonstrates that the SLUG is an 
attractive gain element for the realization of a low-noise 
microwave amplifier. We now consider concrete external 
networks used to embed the device in a 50 fi environ- 
ment. The tasks are to maximize power transfer to and 
from the device and to match the 50 fl source to the 
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FIG. 15. Minimum number of added noise photons in the 
quantum regime riopt versus flux for bias currents h = 1.8 Jo 
(black), 1.9 7o (blue), and 2.0Jo (red). The operating fre- 
quency is 5 GHz and the SLUG parameters are Pl = 1, 
Pc = 0.8, L = 10 pH, and C = 50 fF. 



optimal noise impedance at the desired operating fre- 
quency. For example, to maximize gain we design a con- 
jugate matching network to transform the 50 fi source to 
Ri — jXi. On the other hand, optimal noise performance 
is achieved for an input matching network that trans- 
forms the 50 O generator to the complex optimal source 



impedance Z^^opt = Rs.opt + jXg^opt- Since Rs,opt > Ri 
for typical parameters, it is generally not possible to 
achieve a simultaneous power match and noise match. 
However, it is possible to find a compromise where there 
is reasonable gain and good noise performance over a rel- 



atively broad bias range. Fig. 16 1 shows a schematic di- 
agram of a SLUG-based microwave amplifier with trans- 
mission line matching sections at the input and output. 
To calculate amplifier gain and noise performance, we 
treat the SLUG as a "black box" with scattering and 
noise parameters derived from the calculations of Sec- 
tions IV- VI (Fig. ^). 

As an example, we show in Fig. 17 the frequency- 
dependent gain, noise temperature T„, and added noise 
quanta n for SLUG amplifiers operated with different 
single-section transmission line input couplers with char- 
acteristic impedance in the range from 1-3 fl. Here we 



have used the full expressions (20) and (271 to calcu- 



late the frequency-dependent noise contribution of the 
amplifier in the thermal and quantum regimes, respec- 
tively. The length of the input coupler provides a bare 
quarter-wave resonance at 6.5 GHz; inductive loading 
by the SLUG pulls the operating frequency down to the 
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FIG. 16. (a) Schematic of SLUG microwave amplifier, (b) 
Circuit for amplifier analysis. 
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FIG. 17. (a) Gain, (b) noise temperature, and (c) added 
noise quanta for a 5 GHz amphfier incorporating a 10 pH 
SLUG element with Pl = I, Pc = 0.8, C = 50 fF and h = 
1.8 Jo . The input matching network is a single transmission 
line section with characteristic impedance as indicated. Gain 
and added noise are evaluated at the frequency where the 
quantum noise contribution of the SLUG is minimum. 



desired value of 5 GHz. We remark that the transmis- 
sion line impedances considered here are readily achieved 
with thin-film microstrip technology: for example, a trace 
width of 10 /im and a dielectric with e^ = 4 and thickness 
100 nm corresponds to a characteristic impedance of 2 fJ. 



In Fig. 18 we consider the frequency-dependent gain 
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FIG. 18. (a) Gain, (b) noise temperature in the thermal 
regime, and (c) added noise in the quantum regime for a 
5 GHz SLUG amphfier. The device parameters are (3l = 1, 
I3c = 0.8, L = 10 pH, C = 50 fF, and h = 1.8 7o. The input 
matching network is a single transmission line section with 
bare quarter-wave resonance at 6.5 GHz and characteristic 
impedance 2 fl. 
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FIG. 19. (a) Gain, (b) noise temperature in the thermal 
regime, and (c) added noise in the quantum regime for broad- 
band amplifiers incorporating a 10 pH SLUG element with 
/3l = 1, /3c = 0.8, h = 1.8/0, and $6 = 0.35 $0. The red 
traces correspond to a three-section input matching network 
with quarter-wave resonances at 5 GHz and with character- 
istic impedances of 24.3 $1, 17.4 fl, and 3.0 Jl, derived from 
numerical minimization of the SLUG quantum noise over the 
band from 4.5 GHz to 5.5 GHz. The blue traces correspond 
to a matching network consisting of three sections with char- 
acteristic impedance 29.8 Q, 7.1 fl, and 1.1 fl followed by a 
series capacitance of 38 pF to tune out the imaginary part of 
the SLUG input impedance at a frequency of 5 GHz. 



and noise performance of SLUG amplifiers operated with 
different fixed single-section input coupling networks. 
Due to the nonvanishing cross spectral density Syj, the 
minimum noise temperature occurs at a frequency that 
is somewhat lower than the frequency of maximum gain. 
For a Zo^i = 2 D, input coupler, we achieve noise within 
50% of the standard quantum limit at a frequency where 
amplifier gain is 15 dB, and noise within a factor of 2 of 
the standard quantum limit at a frequency where gain is 
18 dB. 

Finally, we note that is possible to increase amplifier 
bandwidth significantly by coupling the input signal to 
the device via a multisection transformer with stepped 
characteristic impedances. As an example, we show in 
Fig. 19 the frequency-dependent gain and added noise for 



amplifiers operated with different three-section matching 
networks. Here the length of each transmission line sec- 
tion is chosen to provide a bare quarter-wave resonance 
at 5 GHz, and the characteristic impedances were deter- 



mined by numerical minimization of the quantum noise 
contribution of the SLUG in the frequency range from 
4.5 to 5.5 GHz. 



VIII. DYNAMIC RANGE 

The strong nonlinearity of the SLUG leads to gain 
compression and harmonic generation when the device 
is driven with a large-amplitude signal. It is important 
to verify that the SLUG dynamic range will be sufficient 
for the desired application. In Fig. [20^ we plot nor- 
malized SLUG gain versus signal power coupled to the 
device input over a range of bias parameters for /Sl — 1, 
Pc = 0.8, L = 10 pH and C = 50 fF. These plots were 



generated by solving the SLUG equations of motion ( 10 ) 
with a sinusoidal driving term of varying amplitude. De- 
pending on bias point, the 1 dB compression point occurs 
somewhere in the range from -110 dBm to - 90 dBm, cor- 
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responding to input powers from 10 fW to 1 pW. These 
1 dB compression points are comparable to those seen in 
other SQUID-based microwave ampHfiers [7] and 1-2 or- 
ders of magnitude higher than those achieved with typical 
Josephson parametric amplifiers J^ . Amplifier dynamic 
range is determined by dividing the signal power at 1 dB 
compression by the noise power contributed by the SLUG 
over a given bandwidth. In Fig. [20| 3 we plot SLUG dy- 
namic range; here we have used the zero-temperature 
quantum spectral density for the shunt resistors of the 
SLUG. We find a typical value of 130 dB Hz, correspond- 
ing to a dynamic range of 40 dB in an amplifier band- 
width of 1 GHz. For applications related to dispersive 
readout of qubits in a circuit QED architecture, where 
the focus is on measurement of signals at the level of sin- 
gle microwave quanta in bandwidths of order 100 MHz 
to 1 GHz, the dynamic range of the SLUG amplifier is 
more than adequate. 



IX. EFFECT OF INPUT CIRCUIT ADMITTANCE 

In the above analysis, we have solved for the behav- 
ior of the isolated SLUG element, and then treated the 
device as a "black box" with known scattering parame- 
ters for the purpose of designing appropriate matching 
networks. In reality, the nonvanishing admittance at the 
device input and output will modify the device charac- 
teristics, and a complete treatment must take loading by 
the external circuit into account. The scattering param- 
eters will now depend on the particular realization of the 
matching network, and a full exploration of the space 
of design parameters becomes tedious. However, we find 
that the performance of the SLUG amplifier is not greatly 
affected by the nonvanishing input circuit admittance, 
particularly once modest steps are taken to decouple the 
SLUG element from the higher-order modes of the reso- 
nant input matching network. 

To take into account the admittance of the resonant in- 
put matching network, we modify the junction equations 



of motion ( 10 1 to include an additional term represent- 
ing the current drawn by the input circuit. The circuit 



model is shown in Fig. 21 1. The input transmission line 
of impedance Zq can be exactly modeled as a pair of cou- 
pled, time dependent voltage sources E]^ and Eg. These 
are related to the voltages Vl.s and currents /l,s at the 
two ends of the transmission line as follows: 
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(32) 

where t^ is the propagation delay along the transmis- 
sion line. The input current is then determined by the 
additional differential equation 
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FIG. 20. (a) Normalized gain verus input power for a SLUG 
element with Pl = 1, /3c = 0.8, L = 10 pH, C = 50 fP, and 
h ~ 1.8 /q. The different traces correspond to various flux 
bias points, (b) SLUG dynamic range versus flux for various 
current bias points; the device parameters are as in (a), and 
we assume a zero-temperature quantum spectral density for 
the SLUG shunt resistors. 



Using the modified equations of motion for the junction 
phases, we calculate the dc characteristics of the SLUG. 
The LV and V-$ curves of a 10 pH, /3l = 1 SLUG with 
a 10 GHz quarter-wave input transformer are shown in 
Fig. 



22 i-b. We observe sharp Shapiro step-like structure 



(33) 



at voltages corresponding to Josephson frequencies that 
are integer multiples of the half-wave resonance of the 
input circuit. While quantum fluctuations of the SLUG 
shunts smooth out this structure somewhat, it is clearly 
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FIG. 21. (a) Model for circuit analysis with finite input circuit 
admittance, (b) Amplifier circuit with filter inductor Lf to 
decouple SLUG from modes of the input circuit. 
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FIG. 22. (a) 1-V curves of a SLUG operated with a transmis- 
sion hue input circuit with characteristic impedance Zo = 2 fl 
and bare quarter-wave resonance at 10 GHz for $5 = $o 
(black), 0.25 $0 (blue), and 0.5 $0 (red), (b) V-$ curves of 
the same SLUG for h = 1.8 7o (black), 1.9 7o (blue), and 2.0 Jo 
(red), (c)-(d) As in (a)-(b), respectively, for a circuit incor- 
porating a 60 pH filter inductor Lf to decouple the modes of 
the SLUG from the modes of the input circuit. 



desirable to decouple to the SLUG from the higher-order 
standing wave modes of the input circuit, as these modes 
will limit amplifier dynamic range and lead to excess 
noise. 

To suppress the resonances of the input circuit, we 
insert a filter inductor Lf of order tens of pH between the 



FIG. 23. (a) Gain and (b) added noise in the quantum regime 
for SLUG amplifiers calculated using the "black box" scatter- 
ing parameters of the isolated SLUG (red traces) or by solving 
the full circuit model of Fig. |21[ The SLUG parameters are 
/3l = 1, /3c = 0.8, L = 10 pH, C = 50 fF, /;, = 1.8 Jo, and 
^b ~ 0.35$o. The input matching network consists of a 2 f2 
transmission line section with bare quarter-wave resonance at 
10 GHz followed by a filter inductor Lf — 60 pH. 



input coupler and the SLUG element, as shown in Fig. 
[2H 3. In Fig. [22j;-d we plot the SLUG characteristics with 
a 60 pH filter inductor in place. We see that the resonant 
structure is greatly suppressed. 

We can now calculate the gain and noise properties 
of the complete circuit of Fig. |21[ 3 by performing a full 
integration of the amplifier equations of motion. Power 
gain and bandwidth are determined by driving the am- 
plifier with a sinusoidal input tone and monitoring the 
SLUG output at the excitation frequency. In Fig. [23^ 
we plot frequency-dependent gain for the SLUG circuit. 
The blue trace is the result of the full circuit simulation, 
where we have taken a transmission line input with char- 
acteristic impedance Zq = 2 and a length corresponding 
to a bare quarter-wave resonance at 10 GHz, significantly 
higher than the amplifier operating frequency of 4.5 GHz 
in order to compensate for the additional reactive load- 
ing by the filter inductor. The red trace was obtained 
by treating the SLUG as a "black box" with scattering 
parameters calculated as described above in Section IV. 
The agreement with the full circuit simulation is good, 
confirming that the filter inductance has effectively iso- 
lated the modes of the SLUG and the input circuit. 

To calculate the frequency-dependent noise tempera- 
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ture Tn{f), we simulate a "hot load / cold load" exper- 
iment where we compare the power spectra Sv,coid and 
Sv,hot at the device output for source resistances at tem- 
peratures T — and T^, respectively. In the thermal 
regime, we have 



Tnif) 



Sv,coldif) 

Svjiotif) — Sv,cold{f) 



n. 



In the quantum regime, we find 

C0th[hf/2kBin + Tr,)] ^ Sy.hot 

coth {hf/2kBTn) Sv,coid ' 



(34) 



(35) 



The added noise number is then obtained from Eq. |26| 
In Fig. |23Jd we plot the added noise of a 5 GHz SLUG am- 
plifier calculated with the full circuit model and with the 
"black box" scattering parameters of the isolated SLUG. 
The noise magnitude is similar in the two cases, although 
the full circuit solution predicts a higher frequency for the 
minimum in the amplifier noise contribution. We under- 
stand the shift in the frequency-dependent noise charac- 
teristics to be due to a modification of the circulating 
current spectral density Sj by the nonvanishing admit- 
tance of the input network. 



X. HOT ELECTRON EFFECTS 

At millikelvin temperatures electrons decouple from 
the phonons, and the electron temperature of the SLUG 
shunts can be significantly higher than the bath temper- 
ature. Wellstood et al. showed that the electron temper- 
ature Te in a metal thin film resistor is given by 



T, = {p/j:n + T^) 



5U/5 



(36) 



P is the power deposited in the resistor, 
a materials parameter equal to approximately 



fl is the normal metal volume, and T, 



p 



where 
S is 

2 X 10^ W/m^K^ 
is the phonon temperature [23] . The elevated tempera- 
ture of the shunt resistors translates directly to elevated 
noise temperature of the amplifier. For a device with 
fixed /3c, the power dissipation in the shunts scales as 
1/B?. Hot electron effects will be particularly relevant 
for the microwave amplifiers discussed here, as optimal 
performance is achieved for small SLUG inductance, cor- 
responding to large critical currents and small shunt re- 
sistances. 

A proven strategy to promote thermalization of the 
SLUG shunts at millikelvin temperatures is to fabricate 
large-volume normal metal cooling fins in metallic con- 
tact with the resistor element. At low temperatures, the 
inelastic diffusion length is of order several mm |23j ; the 
cooling fins thus allow hot electrons generated in a local- 
ized region of the shunt resistor to diffuse over a large 
volume and thermalize with cold electrons and phonons. 
Wellstood et al. demonstrated a significant reduction 



in the electron temperature of dc SQUIDs incorporating 
400 X 400 iMT? CuAu cooling fins with thickness around 
1 /^m, with measured electron temperatures under 40 niK 
pi] . A similar approach has been used to suppress hot- 
electron effects and reduce the noise temperature of mi- 
crostrip SQUID amplifiers operated in the radiofrequency 
regime 4_. It will be straightforward to integrate nor- 
mal metal cooling fins with area of order 1 mm^ into a 
standard microwave SLUG amplifier geometry without 
compromising the microwave integrity of the circuit. We 
anticipate that the addition of such cooling fins will make 
it possible to attain electron temperatures under 100 mK 
for the device parameters considered here, corresponding 
to operation far in the quantum regime for frequencies in 
the range from 5-10 GHz. 



XI. CONCLUDING REMARKS 

We have presented a comprehensive theoretical treat- 
ment of the SLUG microwave amplifier. Specific advan- 
tages of this approach over competing approaches to low- 
noise microwave amplification are as follows: 

1. The low-inductance device geometry is compact, 
straightforward to model at microwave frequencies, 
and readily integrated into a microwave transmis- 
sion line. 

2. The device input and output are both rea- 
sonably well-matched to a 50 11 transmission- 
line impedance, facilitating broadband operation. 
Moreover, multisection transmission-line input cou- 
plers provide a clear path to attaining bandwidths 
of order GHz while maintaining excellent gain and 
noise performance. 

3. It is straightforward to decouple the SLUG modes 
from the input modes, allowing separate optimiza- 
tion of the gain element and the input matching 
network. 

4. The dynamic range of the amplifier is large relative 
to that required for qubit readout or circuit QED 
applications. 

5. Due to its extremely small magnetic sensing area, 
the SLUG gain element is extremely robust and 
immune to ambient magnetic field fluctuations. 

We believe that we have identified the major techni- 
cal obstacles and outlined a clear path to device opti- 
mization. We anticipate that these amplifiers will be 
attractive in the context of qubit readout in a circuit 
QED architecture [25], either as a near quantum-limited 
first-stage amplifier or as an ultralow noise postamplifier 
following a Josephson paramp. Other possible applica- 
tions include fundamental studies of microwave photon 
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counting statistics [55], or ultralow noise amplification 
for dark-matter axion detection [27] . 
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